Spectral and thermodynamic properties of a strong-leg quantum spin ladder 



D. Schmidiger/ P. Bomllot,^ S. Miihlbauer/ S. Gvasaliya/ C. Kollath,^ T. Giamarchi,^ and A. Zheludev^'* 

^Neutron Scattering and Magnetism, Laboratory for Solid State Physics, ETH Zurich, Switzerland. 
'^DPMC-MaNEP, University of Geneva, CH-1211 Geneva, Switzerland. 
^DPT-MaNEP, University of Geneva, GH-1211 Geneva, Switzerland. 
(Dated: December 20, 2011) 

The strong-leg S = 1/2 Heisenberg spin ladder system (C7HioN)2CuBr4 is investigated using 
Density Matrix Renormalization Group (DMRG) calculations, inelastic neutron scattering, and bulk 
magneto-thermodynamic measurements. Measurements showed qualitative differences compared to 
the strong-rung case. A long-lived two-triplon bound state is confirmed to persist across most of 
the Brillouin zone in zero field. In applied fields, in the Tomonaga-Luttinger spin liquid phase, 
elementary excitations are attractive, rather than repulsive. In the presence of weak inter-ladder 
interactions, the strong- leg system is considerably more prone to 3-dimensional ordering. 



In quantum magnets, the interplay between exchange 
and quantum fluctuations leads to a host of novel phases, 
much richer than their classical counterparts. In par- 
ticular, correlations between the spins can be strongly 
suppressed by quantum effects, leading to quantum spin 
liquid phases with properties quite different from those of 
any conventional ferro- or antiferromagnet. Under mag- 
netic fields these systems undergo quantum phase tran- 
sitions that are akin to Bose-Einstein condensation [1]. 
Among the spin liquids, antiferromagnetic (AF) Heisen- 
berg = 1/2 ladders are the simplest, yet perhaps the 
most important and extensively studied [2]. They com- 
bine the essence of quantum magnetism with peculiar 
features that stem from their one dimensional nature 
[3]. As a result, in applied fields they demonstrate a 
variety of scaling properties, characteristic of the physics 
of one dimensional interacting quantum particles, the so 
called Tomonaga Luttinger liquids (TLL) . Understanding 
which key parameters of the actual spin Hamiltonian con- 
trol these universal features is a formidable challenge that 
requires novel experimental and theoretical approaches. 

In recent years, a general theory of weakly coupled lad- 
ders under strong magnetic fields has emerged [4]. Con- 
siderable experimental progress in understanding strong- 
rung spin ladders was made through the study of the 
compounds IPA-CuCla [5-7] and BPCB [8-13]. Partic- 
ular attention was given to the field-induced quantum 
phase transitions [6-8, 10, 11], and the properties of the 
gapless TLL critical phase at intermediate fields [10, 12]. 

In the case of the strong rung ladder, the spin gap 
in the absence of a magnetic field is already present on 
each rung, protecting the spin-liquid state from the leg 
exchange. A more subtle limit is provided by the strong 
leg (or weak rung) ladder. In that case the existence 
of a spin liquid state is far from obvious, and results [2] 
from an Haldane gap mechanism [14]. This leads to some 
similarities between the two limits but of course also to 
important differences, in terms of the origin of the spin 
gap, excitation spectrum, and the TLL mapping. On 
the experimental side, this interesting problem remained 
elusive since only few studies are available. 



In this paper we report both experimental and theo- 
retical studies of the prototypical strong-leg spin ladder 
material DIMPY [15, 16]. We determine both its thermo- 
dynamic properties and the neutron scattering spectrum, 
and show how to use these data to determine the TLL 
parameters in the strong-leg case. A remarkable quan- 
titative agreement with DMRG calculations gives us a 
precise description of the material, needed to understand 
long range ordering detected at low temperatures. 

The magnetic properties of DIMPY originate from lad- 
ders formed by S = 1/2 Cu^+ ions, that run along the a 
axis of the monoclinic crystal structure [18]. We model 
this compound by the AF Heisenberg two-leg spin ladder 
Hamiltonian 

'H = Jleg ^ Sij-Si^ij-\-Jrung ^ Si^i- Si^2— 9 f^sH Sfj . 

Here Jieg and Jrung are the couplings along the leg and 
rung, respectively, gii^H is the uniform Zeeman field, 
and Si^j are the spin operators acting on site / of the 
leg j of the ladder. At = 0, the ground state of 
DIMPY is a non-magnetic spin singlet separated from 
the lowest-energy triplet excited states by an energy 
gap of A = 0.36 meV [15, 16]. Previous studies sug- 
gested that the application of a magnetic field at T ^ 
leads to a quantum phase transition to the TLL state 
at Hci = 2.85 T [15]. Inelastic neutron scattering mea- 
surements of the dispersion relation for triplon excita- 
tions yielded an estimate of the ratio of exchange con- 
stants as Jieg/^rung ^ 2.2, through a comparison with 
theoretical results obtained with the PCUT method [15]. 
A more detailed measurement over the whole Brouillon 
zone confirmed that the spin Hamiltonian is symmetric 
with respect to leg permutation [16]. This feature allows 
one to conveniently describe the spin dynamics in terms 
of separate antisymmetric (leg-odd "-") and symmetric 
(leg-even "+") structure factors 

S^^\q,i^) (X Y,\{\\S±{q)\0)\''5{uj + Eo- Ex), 

A 

respectively. Here |0) denotes the ground state of H with 
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FIG. 1. (color online) Dynamic spin structure factor of 
DIMPY. (a) Antisymmetric channel. The false-color plot 
shows the DMRG result^. Symbols are experimental data 
for single-triplon dispersion from Ref. [16]. (b) Symmetric 
channel. False-color plot as above. Symbols are positions of 
peaks in inelastic neutron scattering scans shown in (c)-(e) 
(symbols). The white dashed lines in (b) are the limits of 
the two triplons continuum. The solid red line in (c)-(e) is 
the DMRG result scaled by an arbitrary factor and convolved 
with the resolution function of the 3- axis spectrometer [17]. 
The green arrow in (c)-(e) is the lower edge of the two triplons 
continuum. 

^ The oscillations are numerical artefacts. 
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FIG. 2. (color online) Magnetization induced in DIMPY as a 
function of applied field. Symbols are the experimental data 
obtained for if||a at T = 500 mK. The solid line is a T = 
DMRG calculation. The inset shows a blowup of the low-field 
region. 



energy Eo, S±{q) = ^-'^^^(^^,1 ± 5z,2), a the lattice 
constant, and is the sum over all eigenstates |A) of H 
with energy Ex. The two channels can be independently 
probed by inelastic neutron scattering experiments. 

To validate the spin Hamiltonian and to obtain a more 
accurate estimate of the Jieg/^mng ratio, we fit the ex- 
perimental results of Ref. [16] for the full single-triplon 
dispersion present in S''~\q,uj) with quasi-exact numer- 
ical results adjusting both Jieg and Jrung- The calcula- 
tions were performed using the time-dependent DMRG 
method [19, 20] (for specific details see [21]). An almost 
perfect agreement with experiment is obtained over the 
whole Brillouin zone (BZ) with Jieg = 1.42(6) meV and 
Jrung = 0.82(2) meV shown in Fig. 1(a) [22]. The ex- 
cellent agreement with data validates that H is a faith- 
ful description of the system, and that additional terms 
(anisotropics, Dzialoshinski-Moryia etc.) if present are 
extremely small. We obtain Jieg/ Jrung = 1.72(6) for 
DIMPY, which is notably less than the value quoted in 
Ref. [15]. The main difference occurs in Jieg, and we at- 
tribute this difference to the approximation within the 
PCUT method [15] and to fitting in the whole BZ in our 
case. 

The obtained exchange constants were used to cal- 
culate the symmetric structure factor S^~^\qjUj). In 
the strong-rung limit these excitations are attributed to 
multi-particle states with an even number of triplons 
[23, 24]. Assuming no interactions between excitations, 
one expects to see a diffuse continuum of two-triplon scat- 
tering with a maximum of the lower boundary at the cen- 
ter of the BZ ka = tt. Interactions lead to two-triplon 
bound states [25]. In the strong rung limit, these only 
exist below the continuum in a narrow range close to 
the BZ-center. The actual calculated symmetric spec- 
trum for DIMPY is shown in Fig. 1(b) and deviates from 
this simplistic picture. Similarly to the isotropic point 
(Jieg = Jrung) [26], the coutinuum has a local minimum 
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FIG. 3. (color online) TLL parameters as a function of field- 
induced magnetization. DMRG results determined following 
Ref. [21, 29] for DIMPY with Jieg/Jrung = 1.72(6) are shown 
as solid lines and the strong-rung coupling limit as dashed 
lines [21]. The symbols in (b) are extracted from the experi- 
mental heat capacity measurement. 

in the center of the BZ, where most of the spectral weight 
is concentrated. A characteristic 'hat' on top of the 
continuum can be identified. However, the most promi- 
nent feature is a long-lived excitation below the bound- 
ary of the continuum, stable across most of the BZ, at 
0.8 • 27r > ka > 0.2 • 27r. Numerically, integrating the 
singular and non-singular parts of the dynamic structure 
factor up to 5 meV, we estimate that 56% of the spectral 
weight is contained in single-triplon excitations and 14% 
in two-triplon bound states. 

The theoretical results were tested in inelastic neutron 
experiments at the TASP 3-axis spectrometer at Paul 
Scherrer Institute, using the same deuterated single crys- 
tal samples and experimental conditions as in Ref. [16]. 
Typical const ant-g scans measured at T = 1.5 K at sev- 
eral wave vectors that correspond to S''~^\q^uj) are shown 
in Fig. l(c-e) in symbols [27]. The only adjustable pa- 
rameter is an overall scale factor. The quantitative agree- 
ment between theory and experiment is a spectacular val- 
idation of our approach. In particular, it is possible to 
experimentally separate the bound state from the contin- 
uum. This is more delicate in strongly dimerized com- 
pounds [28], or those with large energy scales [24]. 

The fitted Hamiltonian also allows us to interpret bulk 
magnetometric experiments. The measured magnetiza- 
tion curve [30], for a field applied along the a axis at 
T = 500 mK is in excellent agreement with DMRG re- 
sults as shown in Fig. 2. The small discrepancy in the 
very vicinity of Hc^ being due to finite-T effects. The 
onset of magnetization signals the gapless TLL regime. 



FIG. 4. (color online) Inset: Magnetic specific heat measured 
in DIMPY in fields applied along the b axis. Main panel: 
Field-temperature phase diagram of DIMPY. The area be- 
tween Hci and Hc2 is the ordered state. Circles are positions 
of lambda- anomalies in specific heat. The solid line is the 
DMRG result using the adjustable parameter nJ^F- 



Here, the low- frequency long wavelength correlation func- 
tions and other properties are expected to have a univer- 
sal form determined by the so-called Luttinger parameter 
which defines the powers of the algebraic correlations, 
and the velocity u of the linear excitation spectrum. 

In DIMPY, these field dependencies are markedly dif- 
ferent from those in the strong-coupling limit as shown in 
Figs. 3(a,b). In the latter case, K decreases beyond i^ci, 
and returns to unity at saturation Hc2- Throughout the 
TLL phase K < 1^ and the elementary spin excitations 
(spinous) are repulsive. Not so in the strong-leg ladder. 
For DIMPY we see that K increases beyond Hd and 
remains greater than unity at higher fields. This signi- 
fies an attractive interaction between spinous [4]. In the 
direct proximity of saturation at Hc2^ ~ 1, which cor- 
responds to non-interacting spinous. The velocity u in 
DIMPY also behaves quite differently compared to the 
strong-rung coupling case, showing a strongly asymmet- 
ric behavior. This behavior of the velocity will have a 
strong influence on numerous quantities defined by low- 
energy excitations, such as the low energy continuum in 
the gapless phase. As a consistency check, we estimated 
the velocity additionally from the specific heat measure- 
ments discussed below, using the relation C{T) = ^g^^ , 
where C is normalized per spin [31]. This estimate (sym- 
bols in Fig. 3b) are in good agreement with our calculated 
velocity, in particular, considering that the determination 
by the specific heat can be inaccurate, as detailed in [21]. 

TLL physics is endemic to one dimension. Ironically, 
one of the most accurate ways to probe its properties is to 
study the quasi-lD case of weakly coupled ladders. Inter- 
ladder interactions result in three-dimensional long range 
AF ordering at a finite temperature. Assuming unfrus- 
trated and weak couplings, the problem can be treated 
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in the framework of the chain-mean field (MF) theory 
[32] . The characteristics of the ordered state are entirely 
defined by the TLL properties of isolated ladders, with 
only one added parameter: the effective inter-ladder cou- 
pling constant nJ^p. (The form suggests equal coupling 
strength J^p to n ladders). In particular, the field depen- 
dencies of the ordering temperature is given by Eq. (2) 
of Ref. [10]. In this formula the quantity is the am- 
plitude of AF correlations in isolated chains (Fig. 3(c)). 
Combined with the field dependence of (S'^), this gives 
us the field-temperature phase boundary shown in a solid 
line Fig. 4. 

DIMPY was previously hailed as an almost perfect ID 
system that even oX H > Hd remains disordered [15]. 
In fact, more careful specific heat measurements reveal 
a weak but well defined lambda anomaly that appears 
for H > Hci. This can be interpreted as the onset of 
3D long-range order. Typical spin specific heat data col- 
lected in protonated samples for H\\a^ are shown in Fig. 4 
[33]. At each field, the putative ordering temperature Tc 
was identified with the peak position. It is plotted against 
field in symbols in Fig. 4 (right axis). The experimen- 
tally measured phase boundary is in excellent agreement 
with the chain-MF prediction assuming an un-frustrated 
inter-ladder coupling of n J^p = 6.3 /ieV. This agreement 
lends credence that the singularity seen in specific heat 
is indeed associated with the 3D ordering. 

As MF neglects the quantum fluctuations between the 
ladders, J^p may underestimate the real coupling J^ 
This said, given almost the same inter-ladder MF cou- 
pling as in the strong-rung material BPCB (nJ^p = 
6.9 /ieV [10]), the ordering temperature is considerably 
enhanced in the strong-leg case of DIMPY. This effect 
is principally due to the rapid growth of transverse cor- 
relations as defined by Ax^ and their slow falloff due to 
the large showing again the differences between the 
strong leg and strong rung limits. 

This work is partially supported by the Swiss Na- 
tional Fund under MaNEP and Division II. We thank T. 
Yankova for her involvement in the synthesis of DIMPY 
samples. 

Note added: During the final stage of this work 
we became aware of the study by Ninios et al, 
arXiv:1110.5653vl, which contains experimental data 
similar to those shown in Fig. 4. 
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